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Abstract —This paper proposes a data-driven approach to de¬ 
tect the switching actions and topology transitions in distribution 
networks. It is based on the real time analysis of time-series 
voltages measurements. The analysis approach draws on data 
from high-precision phasor measurement units (/rPMUs or syn- 
chrophasors) for distribution networks. The key fact is that time- 
series measurement data taken from the distribution network has 
specllic patterns representing state transitions such as topology 
changes. The proposed algorithm is based on comparison of 
actual voltage measurements with a Ubrary of signatures derived 
from the possible topologies simulation. The IEEE 33-bus model 
is used for the algorithm validation. 

I. Introduction 

Different tools have been developed and implemented to 
monitor distribution network behavior with more detailed and 
temporal information, such as SCADA, smart meters and line 
sensors. Creating observability out of disjointed data streams 
still remains a challenge, though ||Tj. The cost for monitoring 
systems in distribution networks still remains a barrier to 
equipping all nodes with measurement devices. To some ex¬ 
tent, a capable distribution state estimator can compensate for 
the lack of measurement data to support system observability. 
However, topology errors will easily downgrade state estimator 
accuracy. Topology detection is a key component for different 
real-time operation and control functions. Most of literature 
on topology detection is based on state estimator (SE) results 
and measurement matching with different topologies. In Q 
authors propose a state estimation algorithm that incorporates 
switching device status as additional state variables. A nor¬ 
malized residual test is used to identify the best estimate of 
the topology. SE-based algorithms are easy to implement, but 
their accuracy is limited to that of the state estimator. They 
are also sensitive to measurement device placement. In 0. 
the authors provide a tool for choosing sensor placement for 
topology detection. Given a particular placement of sensors, 
the tool reveals the confidence level at which the status of 
switching devices can be detected. Authors in pj are focused 
on estimating the impedance at the feeder level. However, even 
a perfect identification of network impedance cannot always 
guarantee the correct topology, since multiple topologies could 
present very similar impedances. 

In this paper, a real time topology detection algorithm is 
proposed based on time series analysis of phasor measurement 
unit (PMU) data. This approach is inspired by high-precision 
phasor measurement units for distribution systems, called 
micro-synchrophasors or (/i-PMU), with whose development 
the authors are involved 0 - The main idea derives from the 


fact that time-series data from a dynamic system show specific 
patterns regarding system state transitions, a signature is left 
from each topology change. The algorithm is based on the 
comparison of the trend vector, built from system observations, 
with a library of signatures derived from the possible topology 
transitions. The topology detection results are impacted by 
load uncertainty and measurement device accuracy. Therefore, 
the analysis takes load dynamics and measurement error into 
account. The topology detection accuracy is also depends on 
the number of /r-PMUs. But, the simulations shows topology 
detection is converge robustly even with limited measurement 
devices. 

H. Distribution Network Model and Physical 
Topology 

Given a matrix W, we denote its element-wise complex 
conjugate by W , its transpose by and its conjugate 
transpose by W*. We denote the matrices of the absolute value, 
of the real and of the imaginary part of W by |1U|, SRjkP) and 
by 3(tU) , respectively, and with |1U| . We denote the entry 
of W that belongs to the j-th row and to the fc-th column by 
[kPjjfc. Given a vector v, will denote its j-th entry, while 
\v\-j the subvector of v, in which the j-th entry has been 
eliminated. Given two vectors v and w, we denote by (u, w) 
their inner product v*w. We define the column vector of all 
ones by 1 We associate with the electric grid the directed 
graph Q = {V^£), where V is the set of nodes (the buses), 
with cardinality n and S is the set of edges (the electrical 
lines connecting them), with cardinality w\ the set S of the 
switched deployed in the electrical grid, with cardinality r and 
the set V is the set of the electrical grid nodes endowed with 
voltage phasor measurement units (PMUs), with cardinality p. 
Let A G {0, be the incidence matrix of the graph (/, 

A = [a^ ... a^] where aj is the j-th row of A, whose 

elements are all zeroes except for the entries associated to the 
nodes connected by the j-th edge, for which the elements equal 
-1-1 or —1, respectively. If the graph Q is connected (i.e. for 
every pair of nodes there is a path connecting them), then 1 
is the only vector in the null space ker A, 1 being the column 
vector of all ones. In this study, we limit our study to the steady 
state behavior of the system, when all voltages and currents 
are sinusoidal signals waving at the same frequency wq. Thus, 
they can be expressed via a complex number whose magnitude 
corresponds to the signal root-mean-square value, and whose 
phase corresponds to the phase of the signal with respect to an 
arbitrary global reference. Therefore, x represents the signal 
x(t) = \x\s/2s\Ti{ujQt + Ax). 


We will denote the vector of the voltages as u S C", the 
vector of the currents as z S C", and the vectors of the powers 
as s = p + iq G C”, with p,q G K" are the active and the 
reactive power injected at node v. The state of the switches 
is tr e [0,1]’", where \a\i = 0 if the switch v is open, [cr]; = 
1 if the switch I is closed. The measured grid voltages are 
collected iny G C^. We define the trend vector S{ti,t 2 ) G C^, 
as the difference between phasorial voltages taken at the two 
time instants ti and t 2 - i.e. 6 {ti,t 2 ) = u{ti) — u(t 2 )- We 
assume that the deployed PMUs in the distribution network 
take measurements at the frequency /. 

We consider a topology which switches status are 
described by a. Its bus admittance matrix Y is defined as 


[Y.U 


VL:j^k^3k. if j = fc 
Tjfc, otherwise 


( 1 ) 


where Yjk is admittance of the branch connecting bus j and 
bus fc, we neglect the shunt admittances. From Q we see that 
Y„ is symmetric and it satisfies 


= 0 , 


( 2 ) 


i.e. 1 belongs to the Kernel of Y Furthermore, it can be 
shown that if Q, the graph associated to the electrical grid, is 
connected, then the kernel of Y^ has dimension 1. 


We model the substation as an ideal sinusoidal voltage 
source {slack bus) at the distribution network nominal voltage 
Un, with arbitrary and fixed angle </>. We consider, without loss 
of generality, ^ = 0. We model all nodes except the substation 
as constant power devices, or P-Q buses. The system state 
satisfies the following equations 


i = Ya-U 


(3) 

II 

o 


(4) 

Uyty - Py tQy 

V f 0 

(5) 


The following Lemma |j^ introduces a particular and useful 
pseudo inverse of Y^ for our topology detection algorithm. 


Lemma 1: There exists a unique symmetric, positive 
semidefinite matrix G C"^" such that 


lX<,lo = 0. 


( 6 ) 


Applying Lemma from © and Q we can express voltages 
of the grid as a function of the currents and of the nominal 
voltage 

u = X„i-\-lUN (7) 

The following proposition ( ||^) provides a approximation of 
the relationship between voltages and powers. 

Proposition 1: Consider the physical model described by 
the set of nonlinear equations Q, 0, 0 and 0. Node 
voltages then satisfy 

u = + + (8) 

(the little-o notation means that limt7„_>.oo = 0)- 

Equation 0 is derived from a first order Taylor expansion 
w.r.t. the nominal voltage Un of the equation relates powers 


Load differences 


Table I. 


House 1 
House 2 
House 3 
House 4 
House 5 
Aggregate 


Mean (kW) 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


SD (kW) 
0.045 
0.070 
0.113 
0.110 
0.046 
0.184 


Table II. AGGREGATE LOAD DIFFERENCES FOR DIFFERENTS 
FREQUENCY 



Mean (kW) 

SD (kW) 

/ = 1 Hz 

0.000 

0.184 

/ = 0.2 Hz 

0.000 

0.425 

/ = 0.1 Hz 

0.000 

0.604 


and voltages. The approximate solution of 0 has been already 
used with success in state estimation 0, Volt/Var optimization 
0, and the optimal power flow problem Q. 

There is always some noise associated with PMUs, i.e. the 
output of our PMU placed at bus j is 

yj =Uj + ej (9) 


where Cj € C is the error caused by the measurement device. 
A common index for measurement error is the total vector 
error (TVE) pO] . In this paper we assume that the loads have 
constant power factor, and consequently 


[ Pih)]j ^ [ p { t 2 )]j 
[ 7 (^ 2 )]/ 


( 10 ) 


Eurthermore, loads have dynamic behavior, described by 


p{t + 1 ) = p { t ) + np { t ) (11) 

where np{t) is a Gaussian random variable, np(t) ^ 
JV{0, crpp\j). A load measurement data set for five residential 
houses in the Texas, U. S. has been analyzed to drive the 
statistical load model. Some smart meters can measure loads 
every second or a couples of seconds. Load demand (kW) are 
recorded every seconds for a week. Statistical analysis of load 
variations between two consecutive seconds is presented in 
Table H] In the United States, a number of houses are connected 
to one distribution transformer. Therefore, the aggregated loads 
for five houses are considered as the reference for load 
variability in this paper. Lower measurement sampling time 
leads to higher uncertainty in load data variability. In Table 
[H] the aggregate characterization for different frequencies is 
reported. 


III. Identification of Switching Actions 

The basic idea behind our proposed approach is that 
changes in switching status will create specific signatures in 
the voltage waveform measurements. In order to develop the 
theoretical base for the proposed algorithm and its ease of 
mathematical proof, we make the following assumptions. 

Assumption 1: All the lines have the same resistance over 
reactance ratio. Therefore, Gs{Yjk) = aft{Yjk),yYjk. 

Assumption 2: Only one switch can change its status at 
each time. 

Assumption 3: The graph associated to the electrical net¬ 
work is always connected, i.e. that there are no admissible state 
in which any portion of the grid remains disconnected. 





Assumption 4: The initial switches status are known. 


Assumption [T] will be relaxed in Section VI in order to test 
the algorithm in a more realistic scenario. However, it allows 
us to decompose the bus admittance matrix as follow. 


= UErU* + iU^iU* (12) 


where S/j, S/ are diagonal matrices whose diagonal entries 
are the non-zero eigenvalues of 3('Kcr(t-i)) and 
U is an orthonormal matrix that includes all the associated 
eigenvectors and S/ = aT,R. From (|^, it can be showed that 
U spans the space orthogonal to 1. Furthermore, we have 

= + (13) 

with r = (/ — ICq ). Assumption 0 is reasonable for the 
proposed algorithm framework; it worE on a time scale of sec¬ 
onds, and typically the switches are electro-mechanical devices 
and their actions are not simultaneous. Finally, Assumption]^ 
is always satisfied during the normal operation. 

Assume that at time t — 1 the switches status is described 
by a(t — 1) = (Ti, resulting in the topology with 

bus admittance matrix l^cr(t-i)- Applying Proposition IT] and 
neglecting the infinitesimal term, the voltages can be expressed 
as 

= + (14) 

Un 

At time t the £-th switch, that was previously open, changes 
its status. Let the new status be described by a{t) = 02 , 
associated to the topology is Since we are basically 

adding the edge in which switch £ is placed from the graph 
that represents the grid, we can write 

^cr(t) = ^o-(i-i) + (15) 

where yi is the admittance of the line, and ai is the f-th row 

of the adjacency matrix associated with the Since an is 

orthogonal to 1, there exists bn such that Uhn = an. This allow 
us to write 

= (l+i0L)U(T.R + ifi{Yn)bebJ)U* 

= (1 -f ia)-^AU{YR + iR{Yi)btbJ)-^U*A^ (16) 

The voltages satisfy 

u{t)=X,y^^ + lUN (17) 

Un 

From ( [T3] | and ([T^, the trend vector can be written as 

= (18) 

Un 

where 

^ait-i)a(t) = - U{J:r + ^{Yn)bnbJ)-^U* (19) 

We can observe that when there is a switching action, the 
voltage profile varies in accordance to a specific topology 
transition. Since = [a 2 ]-n, for the ease of notation in the 

following we will write as $[o.(t)]_f. The following 

Proposition shows a characteristic of ^[a(t)]_i that is crucial 
for the development of our topology detection algorithm. 

Proposition 2: For every topology transition from the state 
described by a{t—l) to the one described by a{t) by changing 
the switch £, is a rank one matrix. 


Proof: Exploiting ([T^i, ([T^, using Ken Miller Lemma 
GD with some simple computations, we can write 

= pUE^XbJ ( 20 ) 

with ^ 

^ ^ 1 + Tr{^{Yn)bnbJ^R)' 

It’s trivial to see that is a rank one matrix with the 

non-zero eigenvalue Acr(t)_f = associated with 

the eigenvector g[cr{t)]-t = LEj^^fo^and thus can be written as 


The trend vector 6 {t) shows the relationship between switching 
actions and voltage profile. Thanks to Proposition we can 
write it as 


6 {t, t — 1) 




s 


^9(7it)_e 


from which we see that 


S{t) ocrg[„yy_^. (21) 

Therefore, every specific switching action pattern that ap¬ 
pears on the voltage profile is proportional to the eigenvector 
g[a{t)]-e’ irrespective of other variables such as voltages u and 
loads s that describe the network operating state at the time. 
Thus, <7[cr(t)]_{ can be seen as the particular signature of the 
switch action. This fact is the cornerstone for the topology 
detection algorithm in this paper. 


IV. Topology Detection Algorithm 


Assuming the distribution network physical infrastructure 
and the initial switches status are known, we can construct 
a library C in which we collect all the normalized products 
between (/ — egl^) and the eigenvectors for all possible 
switches action 


^a(i-l) = {9[a(t)Yt ■ l)]-r} 

where 

rgk(t)]-. 


( 22 ) 

(23) 


The next step is comparing the trend vector 5{t,t — 1) with 
the entries in the library to identify which switch changed 
its status. The detection process is stated in Algorithm 
The comparison is made by projecting the normalized actual 
trend vector onto the topology library C„(y_iy The 

projection is performed with the inner product, and it allows 
us to obtain the projection index for each vector in 





(24) 


If C[o.(t)]_(, — 1, it means that 5 is spanned by g[a(t)]_^ and 
then that the switch I changed its status. Because of the 
approximation I®), the projection will never be exactly one. 
Therefore, we will use a heuristic threshold, called min_proj, 
based on numerous simulations to select the right switch. If 
projection be greater than the threshold, the associated switch 
is selected. Based on simulations, the min_proj is setted to 
0.98. If there is no switching action, the trend vector will be 












Algorithm 1 Topology Changes Detection 

Require: At each time t, a{t — 1), min_proj = 0.98 
1 : a(t) ■>r- a{t — 1 ) 

2: each PMU at each node j record voltage phasor measure¬ 
ments yj{t) 

3: the algorithm builds the trend vector 6{t,t — 1) 

4: the algorithm projects S{t,t — 1) in the library C-p 
obtaining the set of values 


C = 




1€ •^1’; 


5: if maxC > min_proj then 
6: a{t) ^ argmaxC 

7: end if 


zero as all the C[cr(t)]_,, and the algorithm will not reveal any 
topology transition. Notice that the projection value is used 
to detect the change time too, differently of what proposed 
in m), where instead we used the norm of a matrix built 
by measurements (the trend matrix). With a slight abuse of 
notation, we will say that the maximizer of C is the switches 
status (j such that [a]i = 1 if [a{t)\-i = 0 

or vice-versa \a\i = 0 if [cr(f)]_£ = 1 and c^a(t)]-g its 
the maximum element in C. We tacitly assumed so far that 
all the buses are endowed with a PMU, but this is not a 
realistic scenario fora distribution network. In presence of few 
measurements device the algorithm works the same way. The 
only difference is that we are allowed to take the few voltage 
measures 

y = IvX^^t)-^ + lUN (25) 

U N 

where Ip G [ 0 , 1 ]^^" is a matrix that select the entries of u 
where a PMU is placed, and P is the set of nodes endowed 
with PMU. The trend vector will become 


Therefore measurement noise and load dynamics yield non¬ 
zero values for the trend vector, even if there has not been any 
switching action. The projection index ( [24l i may have values 
near unity, leading to wrong topology detection. 

When a switching action happens, branches of the net¬ 
work are changed and current flows change respectively, thus 
causing abrupt voltages variations. Therefore it helps to avoid 
topology detection errors caused by load uncertainty to con¬ 
sider a proper threshold min_norm for the trend vector norm. 
Moreover the additive noise can make the projection value of 
the trend vector onto the library considerably lower than one, 
even if a topology change occurred. This fact prompts us to 
use a threshold on the maximum projection value, min_proj, 
over which we consider if the trend vector change is due to 
a topology transition. To increase the accuracy of topology 
detection, the following steps are added to the algorithm. 
We assume the ideal case without load and measurement 
uncertainty with the £-th switch change its status at time U. 
Consider the trend vector 

6{t,t-T) = y{t) - y{t - t). 

For t < ti and t > ti + t the projections of the trend vector 
onto the library are all equal to zero, because 

S{t, t-r) = y{t) -y{t-T)={) 

Instead for U < f < U -f t, the trend vector is 

5(Cf-r) = 

leading to a cluster of algorithm time instant of length r 
(or J seconds), in which the maximum projection coefficient 
will be almost one. A possible solution is thus to consider 
a trend vector built using not two consecutive measures, but 
considering measures separated by r algorithm time istants 

S{t,t-T) = y{t) - y{t - t). 


= y{ti) - y{t 2 ) (26) 


The elements of the library vector and their dimension change 
too. In fact one can easily show, using ( |25l ) and retracing ( [T7] i 
and ( |20l i that ( |23| l becomes 




\\Iv'^9[a(t)].i 


(27) 


Of course, if we have only few PMUs, we have to tackle the 
observability problem, i.e. we have to find a way to place 
the PMUs such that we are able to detect topology changes. 
Therefore, we have to minimize number of PMUs and maintain 
the system observability for topology detection. 


V. Measurements and Loads Uncertainty 


So far, we considered the case in which the measurement 
devices were not affected by noise and loads were static. In 
reality, there is some noise associated with PMUs. If we take 
and ([TT|l into account, the trend vector becomes 


5{tl,t2) — Ip{Xcrl^t^) - Xcr(t2))-^jM + et2 - et2 + 


IvX^Ui 

Un 


ti-i 

^ np{t) - inq{t) (28) 
t — t2 


Assume that a topology change has happened at time t when 
we have a cluster of algorithm time intervals of length r (^ 
seconds). The former observations lead to the Algorithm 1^ 
for topology detection with measurements noise and load 
variation. 

VI. Results, Discussions and Conclusions 

We tested our algorithm for topology detection on the IEEE 
33-bus distribution test feeder fl?) , which is illustrated in the 
Eigure[2 In this testbed, there are five switches (namely Si, 
S 2 , S 3 , S 4 , S 5 ) that can be opened or closed, thus leading 
to the set of 32 possible topologies Ti,..., T 32 . Because 
of the ratio between the number of buses and the number 
of switches, some very similar topologies can occur (for 
example the topology where only is closed and the one 
in which only S 2 is closed). In the IEEE33-bus test case, 
Assumption[2about line impedances does not hold, making the 
test condition more realistic. Each bus of the network represent 
an aggregate of five houses, whose power demand is described 
by the statistical Gaussian model ( |TT] i. We tested the entire 
switch monitoring algorithm, in different situations. Eirstly, we 
consider the scenario in which the PMUs are affected by noise 
and the loads are not time varying, and then we add different 
levels of variation to them (associated with different measures 
frequencies). We assume that the buses are endowed with high 













Algorithm 2 Topology Change Detection with Uncertainty 

Require: At each time t, we are given the variables a(t — 1), 
minimizer(t — 1), length_cluster(t — 1) 

1: a{t) ^ a(t — 1) 

2: each PMU at each node j record voltage phasor measure¬ 
ments yj{t) 

3: the algorithm builds the trend vector 6 {t, t — t ) 

4: if \\S{t,t — r)|| < min_norm then 
5: d{t, t — r) 3— 0 

6: minimizer(t) = 0 

7: length_cluster(t) = 0 

8 : else 

9: the algorithm projects S{t,t — t ) in the particular 

library obtaining the set of values 


C = 




I 9[<j{t)]-e £ •^ 1 ’ ) 


10: if maxC >min_proj then 

11: minimizer(t) = argminC 

12: if minimizer(t) = minimizer(<: — 1) then 

13: length_cluster(t) ^ length_cluster(t — 1) H- 1 

14: if length_cluster(t) = r then 

15: a{t) minimizer(f) 

16: end if 

17: else 

18: length_cluster(t) ^ 1 

19: end if 

20 : end if 

21: end if 



Figure 1. Schematic representation of the IEEE33 buses distribution test case 
with the five switches 


precision devices, the p,PMU 03^ affected by Gaussian noise 
such that TVE < 0.05%. It also complies with the IEEE 
standard C37.118.1-2011 for PMUs pO) . Eurthermore, we vary 
the number and position of PMUs, considering the case in 
which every bus is endowed with a PMU, and the case in 
which we have only 7 PMUs deployed, whose placements have 
been chosen experimentally, after Monte Carlo simulations, as 
the one that minimizes the algorithm errors. Eurther research 
is needed to characterize a less onerous and more effective 
placement strategy. The algorithm has been tested in each 
condition via 10000 Monte Carlo simulations The results are 
reported in Table III and Table We can see that, expected, 
the 33 PMUs scenario provides better performances. However 
the results with 7 PMUs are very close, showing the possibility 


Table III. RESULTS AFTER 10000 RUNS WITH 33 PMUs 


SD [kV] 


non 

wrong 

decision 

total 

perc. of 



detections 

detection 

errors 

errors 

errors (%) 

0 


0 

50 

50 

100 

1.00 

0.184, (/ 

^ 1 Hz) 

0 

64 

67 

131 

1.31 

0.425, (/ 

^ 0.2 Hz) 

17 

131 

152 

300 

3.00 

0.604, (/ 

^0.1 Hz) 

72 

211 

249 

532 

5.32 

Table IV. 

Results after 10000 runs with 7 PMUs 

Relative 


non 

wrong 

decision 

total 

perc. of 

SD (%) 


detections 

detection 

errors 

errors 

errors (%) 

0 


0 

56 

56 

112 

1.12 

0.184, (/ 

^ 1 Hz) 

0 

180 

185 

365 

3.65 

0.425, (/ 

^ 0.2 Hz) 

31 

199 

209 

441 

4.41 

0.604, (/ 

= 0.1 Hz) 

76 

245 

298 

619 

6.19 


of a satisfactory implementation of the algorithm also in a more 
realistic framework with few PMUs. Euture developments 
include a deeper study about PMUs placement, better load 
characterization and further, analytic study of the thresholds 
min_norm and min_proj that yield the best performances of 
the algorithm. 
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